arXiv:1502.03610vl [astro-ph.HE] 12 Feb 2015 


5th Fermi Symposium : Nagoya, Japan : 20-24 Oct, 2014 


1 


A Correlation Between Optical, X-ray, and Gamma-ray Variations in 
Blazar 3C 454.3 

Yutaro Tachibana, Nobuyuki Kawai 

Department of Physics, Tokyo Institute of Technology, 2-12-1, Ohokayama, Tokyo, Japan, 152-8550 
Sean Pike 

Deportment of Physics, Brown University, Providence, Rhode Island 02912, USA 
on behalf of the MAXI Team and MITSuME Team 


We present the light curve data of a remarkable blazer 3C 454.3 (z=0.859) in optical, X-ray, and gamma-ray 
bands. Since January 2008, we have been monitoring this object using the 50 cm MITSuME, a optical telescope, 
and detected several flares including extraordinary and simultaneous flares in the 7 -ray and optical bands in 
November 2010. Additionally, the Monitor of All-sky Image (MAXI) has been observing 3C 454.3 continuously 
since August 2009. Using these data and gamma-ray flux observed with Fermi-LAT, we discuss features and 
correlations of flux variations between the energy bands. 


1. INTRODUCTION 


The flat-spectrum radio-loud quasar 3C 454.3 
(z=0.859; Jackson and Brown, 1991) is a well known 
object as one of the most active and brightest sources 
in the gamma-ray sky. Optical and radio observations 
revealed that they are ascribed to strongly Doppler- 
boosted emission from relativistic electrons in a jet of 
plasma. The jet is oriented within ^2° to the line of 
sight (Jorstad et al. 2005). The ’’small and big blue 
bump” has been detected by Raiteri et al. (2007) in 
the optical to ultraviolet band. The small blue bump 
is probably a mixture of a iron lines, Mg 11 lines, and 
Balmer continuum from the broad line region [Ogle et 
al. 2011]. On the other hand, the big blue bump is in¬ 
terpreted as a signature of the thermal accretion disk 
(eg. Pian et al. 1999). The X-to-y-ray radiation from 
3C 454.3 is commonly understood to be produced 
in the plasma jet through an Inverse-Compton (1C) 
scattering process off the same electrons that yield 
the radio to optical synchrotron emission, whereas, 
the origin of the seed photons for the 1C scattering 
is not clear yet. Namely, (i) same synchrotron pho¬ 
tons (SSC), (ii) photons coming from outside the jet 
(EC), or (hi) both of these would be possible to being 
seed photons. The X-ray range is most complicated 
because it might contain the high-energy tail of the 
synchrotron emission besides 1C scattering component 
(Abdo et al. 2010). One example of the SED of 3C 
454.3 is shown in Figure [1] (quoted from Vercellone et 
al. 2010). 

Perhaps one of the most powerful tool to uncover 
the emission mechanism of the jet is time variabil¬ 
ity. Observing the with multi-wavelength for long pe¬ 
riod, and investigating correlations between different 
energy bands, we can build up a physical picture of 
the jet and its surroundings. 

In this paper, we firstly show the light curve data 
in three energy bands, optical. X-ray, and y-ray ob¬ 



log Frequency [Hz] 

Figure 1: The SED of 3C 454.3 during the period MJD 
= 54673-54693. The thin sold, dotted, dashed, 
dot-dashed, and the triple-dot-dashed lines represent the 
accretion disk component, the SSC component, EC on 
the accretion disk, and EC on the broad line region 
(BLR), respectively (Vercellone et al. 2010). The sum of 
those the individual components is drawn by the thick 
solid line. 


tained by MITSuME and SMARTS^, MAXl/GSC^, 
and Fermi/LAT^ respectively, secondly present some 
results of time domain analysis, and Anally discuss 
about correlations of flux variations between these en¬ 
ergy bands briefly. 


^http: //www.astro.yale.edu/smarts/glast/home.php# 
^http://maxi.riken.jp/top/ 

®http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 
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Figure 2: The 7-years light curves of 3C 454.3 from 
MITSuME and SMARTS (R-band), X-ray (2-4 keV) flux 
from MAXI/GSC, 7 -ray (0.1-300 GeV) flux from 
Fermi/LAT. In the shadowed area on the X-ray light 
curve, contaminating by a nearby X-ray source is 
suspected. 

2. THE DATA AND ANALYSIS 

2.1. Light Curves in Muiti Waveiength 

Figure [5] shows the variation of R-band flux from 
MITSuME (Akeno Observatory, red points) and 
SMARTS (Bonning et al. 2012, black points). X-ray 
(2-4 keV) flux from MAXI/GSC, and 7 -ray (0.1-300 
GeV) daily averaged flux from Fermi/LAT. For con¬ 
verting the R-band magnitude to the flux, we assume 
that the 0 mag in R-band correspond to a flux of 
1.42 X 10“^ergcm“^s“^ (Sasada et al. 2011). In the 
middle panel, the gray shaded region indicates the in¬ 
terval which may be contaminated by the nearby X- 
ray source IM Peg (see the ATel #6296). To see the 
X-ray flux variations easily, we summed the daily flux 
data points until satisfying the following condition at 
each data points: Flux/Error > 5. 

Variations, especially activities like flares, are seen 


to be completely correlated across the three energy 
bands. On the other hand, looking closely to varia¬ 
tions in the light curves in the point of view of flare 
amplitudes, the X-ray band seems to show lower vari¬ 
abilities than the optical band and 7 -ray bands. While 
it is certainly affected by a binning size, it was also 
noticed in Vercellone et al. (2010). Additionally, in 
the optical band and 7 -ray band, flux responses origi¬ 
nated from same events (same times) can be different 
from a flare to another. We will present a detailed 
analysis on this point in section 2.4. 

2.2. Color-Magnitude Variation in Opticai 
bands 

Color-magnitude analysis is useful for investigating 
the property of spectral behavior of the optical varia¬ 
tion. We show the R-I color indices as a function of I- 
band magnitude in FigurelH 3C 454.3 shows flux fluc¬ 
tuation on a very short time scales (a few hours), so we 
used the R and I data taken simultaneously using the 
MITSuME tri-color camera for constructing the color 
magnitude diagram. As pointed out by Villata et al. 
(2006), this diagram shows the redder-when-brighter 
behavior which is interpreted as the contribution of a 
stable ’’blue” component, i.e. the thermal radiation 
from the accretion disk (big blue bump) lying under 
the variable ’’red” component i.e. the synchrotron ra¬ 
diation from the relativistic electrons in the jet. When 



I magnitude mjd 


Figure 3: R-I color vs. I-band magnitude for 3C 454.3. 
Golors indicate the date of the observation in units of 
MJD as shown in the right bar. The 

redder-when-brighter trend and the plateau are shown in 
this digram. 
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the jet emission dominates that of the accretion disk, 
a plateau from R-magnitude ~14 towards the bright¬ 
est end will emerge in the color-magnitude diagram. 
It should be note that the magnitude of plateau have 
changed from ~0.65 mag to ^0.5 mag in a time be¬ 
tween MJD ~55500 and MJD ~56800. This may in¬ 
dicate that some physical parameter of the accretion 
disk or the synchrotron emitting electrons (or both of 
these) are dramatically changed within this term. 

2.3. The Ratio of The Optical and The 
Gamma-ray Flux 

Studies about 3C 454.3 have revealed a close con¬ 
nection between the optical and 7 -ray flux behavior 
not only in high activity states but also in non-flaring 
state (e.g. Bonning et al. 2008; Abdo et al. 2010a). 
Assuming photons in optical bands are produced by 
the synchrotron radiation of the relativistic electrons 
in the jet and 7 -ray photons are produced through 
the IC scattering process by the same electrons, we 
can represent observed flux variations in these energy 
bands as 

Aopt - (1) 

^ ( 2 ) 

t thus the flux ratio of the optical band and the 7 -ray 
is predicted as following: 

^ (3) 

F^ 

where is total number of emitting electron, 5 is 
Doppler factor of an synchrotron emitting plasma, ao 
and ag are the spectral indices in the optical and 7 -ray 
bands, respectively, B is magnetic field in the emit¬ 
ting region, and Uext is the external seed photon field 
(Chatterjee et al. 2012). 

Flux ratios calculated using the R-band flux and 
7 -ray flux pairs observed on same days are shown in 
Figure m From this figure,we can notice a trend that 
the flux ratio Fopt /Fy becomes small as the source be¬ 
comes bright. There seems to be some time lag in this 
tendency (the decrease of the ratio lags the increase of 
the brightness), and this point will be investigated in 
section 2.3.1. If the variation is only due to a change 
in Ne, the ratio Fopt/Fy will be constant because the 
ratio does not depend on Ne, as shown in Eq.(3) . 
In addition, as can be appreciated from Eq.(l) and 
Eq.(2), if the variation arose from a change in B or 
UeyA, the optical and the 7 -ray variations will not be 
correlated. Therefore, this supposition is inconsistent 
with the previous research results, except the so called 
’’ 7 -ray orphan” optical-UV flare (eg. Vercellone et al. 
2011). For these reason, this trend may indicates that 
the Doppler factor 5 plays a big role in the flux fluc¬ 
tuation in the optical and 7 -ray bands. The origin of 
flux variation will also be discussed in section 2.5. 




- 11.2 - 11.0 - 10.8 - 10.6 - 10.4 - 10.2 - 10.0 


Log[FDpt(erg ernes')] 

Figure 4: The R-band light curve in the last 6 years 
observation with the corresponding ratio Fopt/Fy as a 
function of time is shown in top and middle panel. In 
this analysis, we selected only the combination data of 
the R-band and the 7 -ray to calculate the ratio between 
these energy bands, hence the number of data point is 
less than the original data (shown in Figure [U. Note 
that the unit of Fopt and Fy is [erg/cm^/s] and 
[photons/cm^/s], respectively, thus the ratio Fopt/Fy is 
not a dimensionless quantity. The bottom panel is the 
ratio as a function of brightness level of the R-band and 
color indicates the date of the observation in units of 
MJD - 50000 as shown in the right bar. 


2.3.1. Z-transformed Discrete Cross Correlation Function 
analysis 

To investigate the time lag between the flux ra¬ 
tio Fopt/Fy and the brightness of the optical band, 
we employed the Z-transformed Discrete Cross Cor¬ 
relation Function (ZDCF) introduced by Alexander 
(1997). This method can correct several biases of 
the discrete correlation function of Edelson & Krolik 
(1988) by using Fisher’s z-transform and equal popu¬ 
lation binning (see the original paper for more detail). 
The calculated ZDCF between the flux ratio Fopt /Fy 
versus the optical band flux Fopt during MJD = 55400 
- 55550 is shown in Figure [5] In this period the most 
prominent change of the ratio and brightness was seen 
as shown in FigurelH A positive time lag indicate that 
the variation of the flux ratio is delayed with respect 
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Figure 5: The R-band light curve in MJD = 55400-55550 
with the corresponding ratio -Fopt/Tly as a function of 
time is shown in top and middle panel. The bottom 
panel is the ZDCF of the flux ratio Fopt/F^ and R-band 
light curve during that interval. The time delay is 
defined as positive if the flux ratio variations lead that at 
R-band flux. 


to that of the brightness, thus the decrease of the flux 
ratio lags by typically about 20 day with respect to 
the increase of the flux in the optical band. The max¬ 
imum likelihood Icr error estimate calculated for the 
points between time-lags = —40 to 40 days is 22.00 
days. The physical meaning of this value is not 
clear at present. In a qualitative manner, however, we 
find the behavior of the flux ratio has the following 
two features: (i) fluctuations on relatively short time- 
scales (a few days) coincides with flare like activities. 



Figure 6: The decomposed light curves of R-band and 
7 -ray into individual flares are drawn by the dotted 
curve. The gray solid curve corresponds to summed flux 
of the modeling light curves, black open circles dente the 
observed light curve, the gray dotted line shows a 
residual flux, and same colors indicate that these flares 
were yielded by a same event in each panel. 


and (ii) a relatively long time-scale decreasing trend 
in active phase (MJD > 55490) in the middle panel of 
Figure m The 22.00 day lag found above is probably 
related to the latter feature, the long-term decreasing 
trend of the flux ratio associated with the large swings 
of the optical flux. 

For the behavior of the flux ratio, described above 
as (i) and (ii), we can interpret this by a existence 
of parameters which change on long time scales (>50 
days) besides the Doppler factor 6. Namely, if the ex¬ 
ternal seed photon field C/ext increase or if magnetic 
field in a emitting region decreases gradually, the flux 
ratio will turn downward as shown in Figure 0] accord¬ 
ing the changes in these parameters. In addition, it is 
also possible that the baseline of the Doppler factor 
progressively increases as an angle between the line 
of sight and the direction of the jet progressively de¬ 
creases. 

2.4. flare analysis 

The fluctuations in the optical and 7 -ray light 
curves of blazars can be interpreted as a superposition 
of individual flares caused by same kind of events lying 
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on a steady baseline flux (e.g. Abdo et al 2010c, Chat- 
tergee et al. 2012). For the optical band, for example, 
the steady component and the flaring component can 
be interpreted as thermal radiation from the accretion 
disk and synchrotron radiation in the jet, respectively. 
In order to investigate properties of the variable com¬ 
ponent, we decompose the light curves into individual 
flares represented by the following function: 

F(t) = Fc-b Fj) -I-, (4) 

where Fc is an assumed constant flux level underlying 
the flare, Fq is the amplitude of the flare, tp is the 
epoch of the peak of the flare, and Tr and Td are the 
rise and decay time of the flare. Fp was constrained 
to be equal to the lowest value of the flux during the 
focusing interval. In our analysis, we assume = Td 
because the acceleration and cooling time scales of 
the relativistic electrons are expected to be substan¬ 
tially shorter than the light-crossing timescale of the 
synchrotron emitting region. Figure [5] shows the op¬ 
tical and 7 -ray light curves for the period containing 
the brightest y-ray flare in the Fermi observation from 
3C 454.3. The individual flares presented by dashed 
lines of same colors in the R-band and the y-ray light 
curves indicate that these variations were produced by 
the same origin. One can see that almost all flares are 
correlated in this period but flux responses in these en¬ 
ergy bands originated from same events may be very 
different. On the other hand, the duration times of 
these flares are roughly equal. These aspects might 
be representing the origin of the observed flares. As¬ 
suming that the sub-equal duration time of the op¬ 
tical and y-ray flares is caused by a light-travel time 
in a same emitting region and the difference of flux 
response is made by a variance of physical parame¬ 
ters in the emitting region such as 5 or F, each flare 
observed might be yielded by different blobs of denser 
plasma possessing almost same size with substantially 
different parameters. 

2.5. Gamma-ray/Optical Correlation 

As shown in Eq.(l) and (2), the optical flux Fopt 
and the y-ray flux are related to each other through 
several parameters. For example, if the flux variation 
is due to a change in total number of emitting elec¬ 
tron Ne, Fry oc Fopt, and if it is due to a change in 
the Doppler factor S, Fj oc ^ inves¬ 

tigate this relationship, we calculate the slope of the 
distribution in the Log(Fopt) vs. Log(F-y) space. Data 
points in the Figure [7] show the observed y-ray and 
the optical fluxes for the two substantially bright pe¬ 
riods: MJD = 55434 - 55532 (upper panel), and 55680 
- 56910 (bottom panel). Color contours indicate the 
density of data points and the lines denote the regres¬ 
sion relation between the fluxes in those two bands. 



- 11.5 - 11.0 - 10.5 - 10.0 - 9.5 

log [R-band Flux (erg cm V^)] 



- 11.0 - 10.5 - 10.0 

log [R-band Flux (erg cm V^)] 


Figure 7: Log(Fopt) - Log(F.y) dependences durring MJD 
= 55434-55532 and MJD = 56800-56910. Color contours 
indicate the density of data points and the solid line 
denote the numerical relation between the fluxes in those 
two energy bands. 


e.g.. Fry oc F^j, where X is the slope of Log(Fopt) 
vs. Log{Fry). In MJD = 55434 - 55532 and MJD 
= 56800 - 56910, parameter X is calculated as 1.79 
and 1.46, respectively. According to Ackermann et 
al. (2010), in MJD = 55060 - 55160, OIR spectral 
index from SMARTS light curve is 1.55 ± 0.05, and 
the y-ray spectral index is 1.5 ± 0.1. For oq = 1.55 
and ag = 1.5, we can obtain the relation F-y oc 
Thus, if we assume that the spectral index was largely 
unchanged and the y-rays are produced by EC pro¬ 
cesses, this result seems to imply that the variation is 
mainly due to a change in S and support the sugges¬ 
tion in section 2.3. 
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3. SUMMARY AND CONCLUSIONS 

In this paper, we present the long-term light curve 
data of 3C 454.3 in three energy bands i.e. optical 
(R-band), X-ray (2-4 keV), and y-ray (0.1-300 GeV) 
provided by MITSuME and SMARTS, MAXI/GSC, 
and Fermi/LAT, and reveal several time variability 
properties using the optical and the y-ray flux data. 
Our main conclusions are as follows. 

1. Using R-band and I-band data from MITSuME, 
the redder-when-brighter behavior and the ’’sat¬ 
uration” effect from R-magnitude ~14 towards 
the brightest end are confirmed by the color- 
magnitude diagram. These results were pointed 
out by some previous research {e.g. Raiteri et 
al. 2008, Zhai et al. 2011). The remarkable re¬ 
sult in this research is the finding of a sign of 
change in the plateau magnitude. The physical 
parameter of the accretion disk or the relativistic 
electron in the jet might have changed gradually 
and significantly. 

2. We investigate a relation between the optical- 
band flux Eopt and the flux ratio of it to the 
y-ray flux Fopt /Fy. The flux ratio becomes de¬ 
creasing as the source becomes brightening, and 
this might indicate that the Doppler factor of 
an synchrotron emitting plasma is playing the 
most important role in the variations in the op¬ 
tical and y-ray flux from 3C 454.3. More de¬ 
tailed analysis is performed for the period of 
MJD = 55400 to 55550. During this period, 
the flux ratio shows the interesting behaviors, 
namely, (i) steep dropping corresponding with 
the flares of the optical and the y-ray, and (ii) 
slow declining during the highly active phase 
(MJD F 55490). It might be interpreted by an 
existence of a change of physical parameters in 
the jet besides the Doppler factor 6 or increas¬ 
ing of baseline of the Doppler factor gradually 
for some reasons. In addition, the ^ 20 days 
time-lag between the flux ratio and the optical 
flux variation is detected. A Interpretation of 
this phenomenon will be addressed in a future 
paper. 

3. We performed a decomposition of the optical 
and y-ray light curves of 3C 454.3 during its 


large flaring activity (MJD = 55490-55540). Al¬ 
most all flares are correlated in this period, 
though flux responses vary significantly from one 
event to the other. This might indicate that the 
size of plasma blobs which radiate the optical 
and the y-ray photons are almost same while 
those parameters besides its size are diverse. 

4. For the two bright periods of 3C 454.3, we inves¬ 
tigate the numerical relationship between Fopt 
and Fj, and F^ oc F^p™ and Fy oc F^pf® are 
obtained in MJD = 55434 - 55532 and MJD = 
56800 - 56910, respectively. These are roughly 
consistent with Fy oc F^pf"^ which is expected if 
the variation is due to a change in the Doppler 
factor S. 
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